Interacting fermions in self-similar potentials 
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We consider interacting spinless fermions in one dimension embedded in self-similar quasiperiodic 
potentials. We examine generalizations of the Fibonacci potential known as precious mean poten- 
tials. Using a bosonization technique and a renormalization group analysis, we study the low-energy 
physics of the system. We show that it undergoes a metal-insulator transition for any filling factor, 
with a critical interaction that strongly depends on the position of the Fermi level in the Fourier 
spectrum of the potential. For some positions of the Fermi level the metal-insulator transition occurs 
at the non interacting point. The repulsive side is an insulator with a gapped spectrum whereas in 
the attractive side the spectrum is gapless and the properties of the system are described by a Lut- 
tinger liquid. We compute the transport properties and give the characteristic exponents associated 
to the frequency and temperature dependence of the conductivity. 

PACS numbers: 61.44.Br, 71.10.-w, 71.30.+h 



I. INTRODUCTION 

The electronic properties of quasicrystalaj have re- 
vealed the importance of the non crystalline order at the 
atomic level. Indeed, the conductivity a of these metal- 
lic alloys displays a unusual behavior since it increases 
when either temperature or disorder increases. It is also 
surprisingly low compared to that of the metals that com- 
posed them. From a theoretical point of view, the influ- 
ence of quasiperiodicity on the spectral and dynamical 
properties of r4JflC|tijon systems has been the subject of 
many studiesafflij'LF . For independent electrons sys- 
tems, it has been shown that the eigenstates, which are 
neither localized nor extended but critical (algebraic de- 
cay) , are responsible of an anomalous quantum diffusion 
in any dimension. Concerning the nature of the spec- 
trum, it depends on the dimensionality but also exhibits 
specific characteristics of the quasiperiodicity. More pre- 
cisely, in one dimension, the spectrum of quasiperiodic 
systems, such as the Fibonacci or the Harper chain, is 
made up of an infinite number of zero width bands (sin- 
gular continuous) whereas in higher dimensions, it can be 
either absolutely continuous (band- like) , singular contin- 
uous, or any mixture. These features are a direct conse- 
quence of the long-range order present in these structures 
despite the lack of periodicity. This absence of transla- 
tional invariance makes any analytical approach difficult 
and one must often have recourse to numerical diagonal- 
ization, except in a perturbative frameworkou. 

Given the complexity of the independent electron prob- 
lem, the influence of a quasiperiodic modulation on an 
interacting system is very difficult to tackle. Attempts 



to solve this ptoblem have been mostly confined to mean 
field solutionsQ or numerical diagonalizationsBE3Eil. We 
have recently proposedFp .a different route, already jOiscd 
with success for periodicE^lo and disordered systemsEjEJ. 
The main idea of this method is to first solve the periodic 
system in presence of interactions ; this is relatively easy, 
either in the one-dimensiagai,|Ca,se for which technique 
to treat interactions existal3^E£lE£P , or even in higher di- 
mensions through approximate (Fermi liquid) solutions. 
In a second step, we study the effect of a perturbative 
quasiperiodic potential via a renormalization group ap- 
proach. Several types of quasiperiodic potentials can in 
principle be studied by this approach but the most inter- 
esting effects come from quasiperiodic potentials which 
have a non trivial Fourier spectr um . Jndeed other po- 
tentials such as the Harper modelElfcH^ who have only 
a single harmonic in their Fourier spectrum are pertur- 
batively equivalent to periodic systemsE^. We have used 
our RG approach to treat interacting smaless fermions 
in the presence of a Fibonacci potentiaOEil. We have 
shown that the existence of arbitrarily small peaks in 
the Fourier spectrum (opening arbitrarily small gaps at 
first order in perturbation) leads to a vanishing critical 
interaction below which the system is conducting. This 
novel metal-insulator transition (MIT) has very differ- 
ent characteristics from those observed in periodic and 
disordered systems for which a finite attractive interac- 
tion is required. These predictions have been successfully 
confirmed by numerical calculationsc^c3. Similar renor- 
malizat|iai|i techniques have been also used in a variety of 
casesESEjEa'^ . Even if some of these properties are spe- 
cific to one-dimensional potentials, these results should 



provide a first step toward the understanding of higher di- 
mensional interacting system in quasiperiodic structures. 
In the present paper, we extend this study to quasiperi- 
odic potentials that generalize the Fibonacci potential. 
We show thaLjthe critical properties obtained in the Fi- 
bonacci caseE3 are generic of other self-similar systems. 
Our results are in agreement with the recerit-inumerical 
results obtained on precious mean potentialsES. The pa- 
per is organized as follows : in Section II, we present the 
model on the lattice and derive its continuous version 
for any potential using a bosonization technique. We de- 
tail the renormalization group treatment of the bosonized 
model and the computation of the flow equations for the 
coupling constants. In Section III, we recall the results 



for the well-known Mott transition (periodic case) and 
we describe the physics of the disordered case for which 
a different kind of MIT occurs. We then discuss the most 
interesting situation : the quasiperiodic case. We explain 
why the non trivial self-similar Fourier spectrum induces 
a MIT whose characteristics are intermediate between 
the periodic and the disordered potentials. The physi- 
cal consequences are discussed in the Section IV with a 
special emphasis on the transport properties. We also 
discuss the question of the strong coupling regime. Con- 
clusions can be found in Section ^ and some technical 
details are given in the appendices. 



II. DESCRIPTION OF THE MODEL AND 
RENORMALIZATION 

A. The Model 

We consider a system of interacting spinless fermions 
in a one-dimensional lattice of linear size La {a — 1 being 
the lattice spacing) described by the following Hamilto- 
nian : 

i 

V^^nini+i+^^WiUi (1) 

i i 

= Ht+Hv+Hw (2) 

where q (resp. q) denotes the creation (resp. annihila- 
tion) fermion operator, rii = cl Ci represents the fermion 
density on site i. In (|l|), t represents the hopping integral 
between sites and V controls the strength of the interac- 
tion between nearest-neighbor particles. In addition, the 
fermions are embedded in an on-site (diagonal) poten- 
tial W . In the following, we consider three main cate- 
gories for W : (i) a simple periodic potential of the form 
Wi — \cos{Qi)] {ii) a random potential uncorrelated 
from site to site; (Hi) a quasiperiodic potential whose 
study is the aim of this paper. In this latter case, we will 
focus on the general class of precious mean potentials 
described in Appendix H. 



In order to treat the interactions in (||), it is conve- 
nient to write the fermion operato rs .. in . t erm of boson 
ones. This bosonization techniquell3ll3't2r provides a 
good description of the low-energy physics of the Hamil- 
tonian Ho{t,V) = Ht + Hy- For completeness and to 
fix the notations we give a brief summary of this method 
in Appendix O. Within this framework, HQ{t,V) can be 
rewritten : 



Hoit.V) = ^j dx{uK){nIif + (-|) {d,4>f 



(3) 



where (j) and II are conjugate bosonic fields respectively 
related to the long wavelength part of the density and the 
current (see Appendix H) . All interaction effects can be 
absorbed in the so-called Luttinger liquid parameters : u, 
the velocity of the charge excitations, and K which con- 
trols the behaviour of the various correlation functions 
(see below). For V ^ t^ analyti c expressi ons of these 
parameters can be obtained (see (B35-B36)). However, 
the bosonic representation is in fact quite general and 
the expression (0) provides the effective Hamiltonian de- 
scribing the low-energy physics of MM-one-dimensional 
interacting spinless fermionic systemE3'E2l. 

Concerning the couphng to the lattice potential W , 
one has, in the continuum limit (see Appendix H) : 



i?i 



w 



dx W{x) p{x) 
dx W{x) X 



(4) 



--d^(j){x) + — coa{2(j){x) ~ 2kFx) ] . (5) 



The various physical observables can be expressed in 
terms of the boson fields. For example, the correlation 
function of the 2kp part of the density is : 

R{x,T,x',t') = {TrP2kF{x,T) p2kp{x\T')) (6) 

^ /Tre*^'^*-'^''^^e^*^'^^'^'''^'A (7) 

where Tr is the time-ordering operator for the imaginary 
time T. In absence of the perturbation Hw, one has : 



R{x, r, x' , t') 
where r = {x,t). 



r — r' 



2K 



for |r- r'l > a, (8) 



B. Renormalization Group Analyzis 

To study the influence of the potential W, we use a 
standard RG approach by analyzing perturbatively the 
renormalization of the correlation function ([7|) computed 
with the full action of the system. First, note that via a 
redefinition of the bosonic field 6 : 



K 

{x) = (j}{x) 

u 



dyWiy), 



(9) 



the term proportional to dx(t> in (0) can be absorbed in 
the quadratic part of the action which thus writes : 



So = — I dr 



;nf(**' 



#1(4-^ 



Introducing the Fourier components by : 

W{x) = A ^ Ty(Q)e"3^ , 
Q 

the potential part of the action reads : 
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u(2'Ka) 



■Y.W{Q) 



(10) 



(11) 



(12) 



Q 



dj. ^i{24>{x.T)+Q-x) j^ ^i{-24>ix.T)+Q+x) 



where Q^ — Q zt 2kp and g = 27raA. Treating 5*^^ per- 
turbatively and imposing that the asymptotic behavior 
(^ is unchanged when varying the cut-off a leads to the 
renormalization of the parameter K and of the Fourier 
components of the potential W. The procedure is de- 
tailed in Appendix O. The RG equations are given by : 



f - --=-")• 



dyq 
dl 



^ - (2-/0yQ, 



with 



s(0 = J E 2^Q [J(QMO) + J(Q-«(0)] 



(13) 
(14) 

(15) 



where the yq — \aW{Q)/u are the dimensionless Fourier 
components of W and / is the scale factor defined by 
ail) = a(0)e' where a(0) = a is proportional to the 
original lattice spacing a. In (|l5|), J is a function whose 
precise form depends on the cut-off procedure used to 
eliminate the short distance degrees of freedom (see Ap- 
pendix O). Different kind of functions are considered 
below, but one typically has : 



j(g*a(0) 



1 for a{l) < \IQ^ 
otherwise. 



(16) 
(17) 



Physically, this means that the renormalization is equiv- 
alent to an investigation of the low-energy properties in a 
window around 2k p in the reciprocal space whose width 
is proportional to e~'. Thus, the ]ull Fourier landscape 
of the potential determines the scaling of the parameters 
K and j/g. This is summarized in Fig. n^. 



III. CRITICAL PROPERTIES 

In this section, we discuss the relevance of the inter- 
actions for different potentials. In other words, we ana- 
lyze the possibility of phase transition when varying the 
strength of the interactions. 



A. Periodic potentials 

This is the simplest case since the potential involves 
only one harmonic : 

W{x)^\cos.{Qx). (18) 

Inserting the foim (|l^) in ( [l3[]l4| ) leads to the following 
flow equational3 : 



dK 

nr 

dy 
dl 



-i^V[J(Q+a(0)+J(Q"a(0)], (19) 
{2-K)y, (20) 



where y = Xa/u. Two different behaviors have to be 
distinguished depending on whether zLQ coincides with 
2k p or not. 



1. Incommensurate case (Q^ ^ and Q 7^ OJ 

In this case, there always exists a scale I* such that 
i{Q'^a[l*)) = 0, as shown on Fig. |l|. At this scale the 
renomalization of K essentially stops. K thus converges 
towards a fixed point K* ^ and the potential W is irrele- 
vant. The ay s te rn remains a Luttinger liquid with gapless 
excitationsllMJ and the correlation functions decay with 
an effective exponent K* . This can be understood by the 
fact that, in this case, the Fermi level does not lie in the 
gap opened by the periodic potential W . 



2. Comm,ensurate case (Q^ = or 1 



0; 



Suppose for instance that Q ~ 2kp. In that case 
i{Q~a{l)) = 1 at all scales and the renormalization of K 
cannot be stopped, as shown on Fig. ||. The potential is 
commensurate and would, for non interacting electrons 
{K — 1) open a gap at the Fermi level. For interact- 
ing particles, the effect of the potential is given by the 
flow ( p^J20| ) whjphrisj now a simple Beresinskii-Kosterlitz- 
Thouless flowEl'^M Thus, for K > 2, the potential W 
is irrelevant, and the system remains a gapless Luttinger 
liquid. For K <2, the potential is relevant and the sys- 
tem flows to strong coupling regime. The strong coupling 
fixed point is not reachable by the perturbative flow but, 
since the system is described in this case by a simple 
sine-Gordon action, we know, fijQj Ti . oth er methods, the 
physical properties of this phaset3't2lE3. A gap opens in 
the spectrum between the ground state and the first ex- 
cited state. An estimate for this gap can be obtained 
from the RG analysis. If I* denotes the lengthscale at 
which y ~ 1, the gap is given by : 



A c>c e 



(21) 



In the case where y ^ (2 — K) one can neglect the renor- 
malization of K in (pOJ ) and we obtain : 
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(22) 
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FIG. 1; Schematic view of the renomialization procedure. The function J (see text) filters the Fourier components W{q) of 
the potential and only keep those in a narrowing window around the Fermi level (i e. for q = 2kF) as the scale increases. We 
have displayed the windows for two different scales h and h > h- For the periodic case the potential has only one peak. If 
the potential is commensurate with 2kF (a), the peak is inside the windows at all scales and the renormalization never stops. 
If the potential is incommensurate (b) it only acts if the peak remains in the windows (which is the case for h but not at 
the larger scale h)- For the disorder (c) case the potential has component at all q thus the strength of the potential efficient 
in the renormalization is directly proportional to the size of the window. For the quasiperiodic case (d) the potential has a 
complicated peaks structure that gives rise to a novel behavior for the MIT. 



Note that for the non interacting point {K =1), one re- 
covers the linear scaling of A with respect to the strength 
of the potential A expected by the first order perturbation 
theory. This MIT induced by the interaction is known 
as the Mott transition. Let us emphasize that this criti- 
cal value Kc ~ 2, separating a metallic phase {K > Kc) 
from an insulating one {K < Kc), corresponds to attrac- 
tive interactions between fermions. 

The later discussion on the single harmonic case still 
holds for a potential with a finite number of Fourier com- 
ponents. The possibility of an insulating regime is then 
offered when the Fermi level is in one of the gaps opened 
by the various frequencies of the potential. The most 
interesting situation thus arises when the Fourier spec- 
trum of the potential is dense or continuous. Such po- 



tentials can be encountered in several physical systems 
but we focus here on two of them : the random and the 
quasiperiodic potentials. 



B. Disordered potentials 

Let us consider a potential provided by a random vari- 
able with uniform probability whose Fourier components 
satisfy : 



W*{Q)W{Q')^X6qq, 



(23) 



Using the general expressions (^_3|) and (JIJ) one obtains : In Eq. ( |24| ) , the sum actually stands for an integral so 

one has : 
dK _ K^y^ 

IT ~ 2~ 



^[j(Q+a(0)+J(g-a(0)],(24) 



l^'--)"- 



(25) 



^[j(Q+a(0) + J(O-a(0)] 



dQ [J(g+a(0) + J(Q"a(0)] ~ 477 ~ e-K 



1 

«"(7) 



(26) 



This means that the renormalization of K is directly pro- 
portionnal to the window width around 2kF at scale I. 
In the limit of weak disorder the RG equations (E4 — 125| ) 
can be integrated neglecting the renormalization of K : 
y{l) = y(0) e^^--^)'. One then has : 



dl 



(27) 



where C is a constant. One deduces the existence of a 
critical point at Kc = 3/2. Our approach thus aeafitalises 
the RG treatment specific to the disorder casetSEj. For 
K < Kc, the system is insulating for any filling whereas 
it can be metallic for sufficiently attractive interaction. 
Note that for K = 1, the system is an insulator as ex- 
pected for a one-dimensional disordered non interacting 
system. By contrast to the periodic case, this transition 
point is independent of kp . Physically, this can be under- 
stood invoking the proximity of several localized states 
(for K = 1) at arbitrarily short distance that can couple 
to each other if the interaction is attractive enough. It 
is interesting to remark that the critical value in the dis- 
ordered case {Kc = 3/2) is smaller than in the periodic 
case {Kc — 2). This means that the localization induced 
by the disorder is more easily destroyed by attractive in- 
teractions than the one induced by a finite width gap. 
In this context, it is clear that dense Fourier spectrum 
are able to provide rich outstanding situations where the 
position of the Fermi level could, in principle, determine 
the value of a critical point. 



C. Quasiperiodic Potentials 

To analyze such situations, we consider quasiperiodic 
potentials. Several types of potentials can play this role 
but, as already explained, the most interesting ones are 
those that have a non trivial dense Fourier spectrum. 
Here, we focus on quasiperiodic potentials obtained by 
substitution rules which lead to quasiperiodic effects even 
at the perturbative level. Among all possible choices, we 
consider the simplest one, known as precious mean po- 
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FIG. 2: Fourier transform of the 15 approximant of the 
Golden Mean (Fibonacci) potential (610 sites per unit cell). 



tentials, that are given by the following iterative scheme : 

A -> A''B, B ^ A. (28) 

We associate, to each site, a diagonal potential that can 
take two discrete values Wa = +A/2 or Wb = —A/2. 
Let us note that a global shift of the Wi always allows 
to deal with a zero-averaged potential so that we can 
set W{Q) = 0. For fc = 1, one recovers the famous Fi- 
bonacci chain associated to the Golden Mean, for fc = 2 
one has the Silver Mean sequence, etc. We give in the 
Appendix H a brief description of these sequences and a 
detailed calculation of their Fourier tranform. As it can 
be inferred from Fig. || (fc == 1) and Fig. || (fc == 2), the 
Fourier spectrum is dense in [0, 27r] iup-the quasiperiodic 
limit, and has a multifractal structureo. As in the pe- 
riodic case, we have to consider several situations since 
we expect a strong dependence of the physical proper- 
ties with respect to the position of the Fermi level. It 
is clear that if there exists a wave vector Q ~ 2fci;' such 
that W{Q) is large compared to the other Fourier com- 
ponents, the electrons behave as if they were embedded 
in a periodic potential, whereas in the opposite case, i. e. 
the Fermi level lies in a very small gap, one could ex- 



IWI 



4 

ft 

X 



.i.i..lI,.,[i,].I...Ij,.i.[.i.I,i.,I.i,j 






FIG. 3: Fourier transform of the 10*'' approximant of the 
Silver Mean potential (1393 sites per unit cell). 



FIG. 4: Behavior of E when the Fermi level 2kF 
close to a dominant peak (Fibonacci potential). 
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pect another behavior. Of course, the notion of far or 
close from a gap has only a sense once we have specified 
the maximum scale Imax up to which we study the RG 
equations. In our case, since we only consider approxi- 
mant (arbitrarily large), a{lmax) is given by the inverse 
of the typical distance between gaps. Beyond this scale, 
the RG flow is no more sensitive to the precise structure 
of the Fourier spectrum. One can then encounter three 
situations : 



1. W{2kF) IS large 

There is an harmonic of W a.t Q ~ 2kp such that 
W{2kp) is large compared to the other Fourier compo- 
nents. In other words, W{Q) opens a finite gap that 
dominates the low-energy physics up to a scale given by 
{Q — 2kp)~^. The function S that governs the flow of 
K is then completely dominated by this component and 
behaves as e*^'*"^^-'' (see Fig. ^). This actually defines 
the proximity of a "dominant" peak. 

In this case, one recovers a critical value Kc = 2 corre- 
sponding to that obtained in the periodic case, separating 
a metallic phase {K > 2) from a insulating one {K < 2). 



2. W{2kF) IS small 

The Fermi level lies far from the main gaps. An exam- 
ple of such a situation is obtained at half filling {2kp = 
tt). In this case, as shown in Fig. |5| S(/) ~ gi'i-^K-2)i 
for any kind of regulator. 

Of course, this is not strictly speaking an exponential 
decrease and there is, in particular, some oscillations (in 
log-log) that are reminiscent of the multifractality of the 
Fourier spectrum. However, one can fairly approximates 
E{1) by e^(^~^'' and we obtain a critical point at Kc = 
1. This critical value corresponds to the non interacting 



point V — 0. This means that when the Fermi level lies 
in zero width gap (small peaks) the slightest attractive 
interaction allows the system to become metallic whereas 
it is insulating at V — 0. From a spectral point of view, 
this means that an attractive interaction close the gap 
and allows for arbitrarily low-energy excitations. This 
point will be discussed more in Sec. IV B. We would like 



to stress that other positions of the Fermi level {2kp = 
0.5 for Fibonacci for example) also gives Kc = 1, and 
that similar situations occurs for other precious mean 
potentials (see Fig.^). 
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FIG. 5: Behavior of H at half-filling {2kF = tt) for the 
Fibonacci potential. We have displayed the results for three 
different regulating function J : J(a;) — e ^ ('^)) ■^{x) = 
1/(1 -I- x^) (o) and the step function J{x) = 1 if a; < 1, 
otherwise (v)- 












y = 0,31214 + -1 ,9796x R= 0,99987 
- y = -1 ,4951 + -1 ,9798x R= 0,9997 
y = -3,4229 + -1 ,9906x R= 0.99098 




regime K > Kc (where Kc = 2 for periodic potential, 
Kc = 3/2 for disordered potentials and Kc — I for spe- 
cial filling factor in quasiperiodic potential), the system 
flows to weak coupling so the RG can be used at arbitrary 
scales. For the insulating side K < Kc the potential flows 
to strong coupling. One can thus use the RG to obtain 
the properties up to the lengthscale I* such that the po- 
tential at this lengthscale is of order one. The behavior 
beyond this lengthscale (or below the correponding en- 
ergy scale) can not be accessible by the RG and should 
be treated by other (non per turba tive) methods. This 
point will be discussed in Sec. IVB, 
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FIG. 6: Behavior of H at half-filling {2kF = tt) for the 
silver mean potential. We have displayed the result for three 
different regulating function J : 3{x) — e~^ (A), J(a;) — 
1/(1 -I- x^) (o) and the step function J{x) = 1 if a; < 1, 



otherwise (v) 



3. W{2kF) is intermediate 



Finally, there are intermediate situations for which S(^) 
does not have an exponential behavior, even approxi- 
mately. In these cases, it is impossible to simply extract a 
critical exponent and one needs a non perturbative treat- 
ment of the problem to determine a possible transition 
point (see Fig. 0). Note that in this case, one still ob- 
serves the oscillations. 
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FIG. 7: Behavior of ^ when the Fermi level 2kF = 0.5 
(Fibonacci potential) is in an intermediate situation. 



IV. PHYSICAL CONSEQUENCES 

In addition to the phase diagram and the critical prop- 
erties, determined in the previous section, the RG al- 
lows to extract many physical properties. In the metallic 



A. Transport properties 

Both the d. c. and a. c. transport properties of the sys- 
tem can in principle be extracted from the Kubo formula. 
The fermionic current operator can easily be written us- 
ing the bosonic variables and reads : 



j = vf{iI)\'4>r.-'iP\'iPl) 
~ uKIi = dr(t)- 



(29) 
(30) 



Thus the conductivity is simply given by the correlation 
function : 



cr(w) 



2uK 



+ X(w) 



(31) 



x(aj) is the retarded current-current correlation function : 



X{'^) = {j\j)u 



dxj di([j(x,t),j(0,0)])e'"*, 

where j is given by (p3) , and L is the size of the system. 
In the absence of Hw , since the Hamiltonian is quadratic 
in the bosonic variables (see (0)), (p3) is trivially com- 
puted and the conductivity is given by : 



cr(w) ~ uK 



n uj 



(33) 



The system is a perfect conductor and uK plays the 
role of the standard plasma frequency. Computing fully 
( pl| ) in the presence of H^ is of course impossible, but 



one can use an hydrodynamic approximatioi 



lajijl 



us- 



ing the so-called memory function formalismEj which is 
well adapted to the one-dimensional situation. For com- 
pleteness, we sketch the main steps of the method in 
Appendix O. The conductivity thus writes : 



a{uj) 



i2uK 



M{uj) 



(34) 



where the function M is perturbatively given by (see Ap- 
pendix O) : 



M{uj) 



[{F-F)l-{F-F)l^,]/u: 

-x(o) 



(35) 
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The F = [j, H] operator takes into account the fact 
that the current is not a conserved quantity and {F; F)^ 
stands for the retarded correlation function of the oper- 
ator F at frequency uj computed in the absence of the 
scattering potential {Hw — 0). Using (||) leads to : 

2m 

F(x) = — W(x) sin(2(j>(x) - 2kFx). (36) 

na 

The memory function can be easi hf com puted for periodic 
potential with a single harmonicE^Eil or for an uncorre- 
lated disordered potential. In order to get the behavior 
beyond the simple perturbati(m.|it| is necessary to couple 
(||) with the RG calculationlla-y. The RG is iterated 
until the frequency or temperature is comparable to the 
renormalized cut-off. The memory function (Bq) with the 
renormalized parameter gives then the conductivity. Us- 
ing the expression of the current pq) one can compute 
the scaling dimension of M. From (|35|) one gets : 



M{uj) = - [ dx f dr(e*"^-l) 



J2\WiQne^^ 
Q 



(7) -^''^ 



Eq. ( |37| ) is essentially the expression of S appearing in 
the RG calculation (p^). The differences are just trivial 
scaling factors. Thus if the function S varies as: 



S(0 oc e(4-2i^-M)'^ 



(38) 



where /i is a real number, then from ( pT\ ) one obtains : 

M{1) ex e^3-2K-^.)t^ ^gg^ 

The function S thus gives directly the conductivity. Stop- 
ping the RG flow when a{l*) ^ max(a;, T), and replacing 
M{1*) in ( |34| ) gives both the d. c. and a. c. conductivity. 
Simplified expressions can be obtained for very weak 
potential, and high enough temperatures or frequencies. 

I 



In that case, one can neglect the renormalization of K 
in the RG flow. Since M is small in this regime, the 
conductivity is given by : 



a(aj, T = 


^0)^1-^^-^-°), (40) 


a{uj = 


T) - (^^) 


^' M(a; = 0,T)' ^^^^ 


leading to : 




aiiu) oc a;2if-5 


a{T) ex T^-"^^ (commensurate), 


a{uj) ex a;2i^-4 


f7(T) ex T^-^^ (disordered). 


aiiu) ex c^2if-3 


(7(T) ex T^-^^ (quasi, for K^ ^ 1) 



A sketch of the conductivity is shown in Fig. g. 

In the metallic regime K > Kc there is, in addition, a 
Drude peak, whose weight is given by the renormalized 
Luttinger liquid parameters as D = u*K*. In the insu- 
lating regime, since the RG can only be pushed until the 
scale Ic for which S(Zc) ~ 1 these expressions, and their 
generalizatiop-jftthen the renormalization of K is taken 
into accountlij^til, are only valid for temperatures and 
frequencies smaller than ae " . This corresponds to the 
Mott gap for the commensurate system and the pinning 
frequency (inverse localization length) for the disordered 
one as will be discussed in more details in Sec. IV B| . 



For a mesoscopic system of size L, it is often more in- 
teresting to compute the conductance G of the system as 
a function of the size L. This is in general a much more 
complicated calculation, specially for interacting systems 
for which is is difficult to use the Landauer formula. How- 
ever, it is possible to extrac t i th p scaling behavior from 
the memory function as wellEZrt^. Indeed, using / — GV 
{V being the electrostatic potential) the conductance is 
given by : 



G= lim 



1 



i^^o i(a; + iS) 



L 



^^^'_dx'y'diV(-+^^)(*-*')(j(x = o,ty,jix',t'))] I 



(42) 



r 



which, rewritten in term of the conductivity, gives : 
G = 



1 /-^/^ 
dq- I dxe*«^a(tj -> 0, q) (43) 



^ J-L/2 

, sm(qL/2) , 



restricting the integral over q. The conductance (in units 
of e^ /h) is thus simply given by : 



G: 



l/L 



dqa{u> — > 0,q). 



(45) 



l/L 



(44) 



When Hw — the conductance is simply G = K. 
Note that the fact that conductance is not G = 1 for 
the pure system but depends on the interactions is in 
The cardinal sine can be considered as a simple cut-off, general an artefact, that has to be corrected depending 



a(T) 



-2K-H) 



FIG. 8: Schematic dependence of the temperature depen- 
dence of the conductivity. Above the MIT point K > Kc the 
conductivity (dashed line) is power-law like with an exponent 
depending on the interactions. On the insulating side this 
behavior remains valid in the weak coupling regime, i. e. for 
temperatures above the crossover scale T^ (see text). 



on the systemua. However, since we focus here on the 
effect of a scattering potential, this is not important for 
our purpose. 



A simple way to compute the conductance in (45) is 
again to iterate the RG flow until the cut-off is of the 
order of the size of the system. Up to that point the q 
dependence can be neglected and cr(w) ^ 1/M{uj). The 
corrections to the conductance are thus given by : 



AGoc - 



1 



LM{l = \n{L/a))' 



Thus, with the scaling (|39|), one has : 
AG oc L2^+^-^ 



(46) 



(47) 



For the disordered case, in the absence of any renor- 
malization of the disorder, (E^) leads to a variation of 
the resistance AR — 1/AG oc L in the absence of inter- 
actions which is nothing but Ohm's law. Interactions and 
renormalizatiosi of the parameter K by disorder change 
this scalingl3E3. Of course the renormalization of dis- 
order also affects the exponents through the renormal- 
ization equations (0-|l3) . The faster increase of the re- 
sistance with size is the sign of Anderson localization. 
As discussed before, the RG results can only be used if 
the renormalized scattering potential remains small com- 
pared to 1. This means that AG < e^/ft,. To go beyond, 
one needs to know the strong coupling fixed point. The 
lengthscale at which this happens is of course the Mott 
length for the commensurate potential and the localiza- 
tion length for the disorder as we now discuss in more 
details. 



the insulating side, on the other hand, the RG flows to 
strong coupling. It thus defines a scale Ic for which the 
function 2 becomes of order one. Below the lengthscale 
L ~ ae'= the RG still gives directly the physical prop- 
erties, as shown for example on Fig H. Note that this 
lengthscale can be quite large if the system is close to the 
MIT or for instance when the Fermi level is far from one 
of the main peaks. Beyond this lengthscale a knowledge 
of the strong coupling fixed point is in principle necessary 
to describe the physics of the system. Fortunately many 
of the properties can still be inferred directly from the 
Hamiltonian. Let us examine the various cases. 

For the commensurate case, we know that the potential 
opens a gap in the spectrum. Thus, we can relate the 
crossover scale to the gap by : 



A/W ■ 



-ic 



f 



(48) 



Thus the RG gives directly the gap. The crossover scale 
is in this case the so-called Mott length ^ above which 
the correlation functions have exponential decay. Sim- 
ilarly the conductivity decreases exponentially for tem- 
peratures below the gap. 

In the disordered case, the situation is more subtle. 
The disorder does not open a gap, but we know that 
if it is strong enough, wavefunctions are exponentially 
localized, with a localization length of the order of the 
lattice spacing. One can thus again relate the crossover 
length to the localization length of the system byllj : 



^ oc ae ' 



(49) 



Because of the exponential localization of the wavefunc- 
tions, £, also defines a scale below which most of the in- 
teraction effects stop to be important. 



B. Beyond the RG 



In the metallic phase the RG gives the full physics 
of the system, up to arbitrarily low energy scales. On 
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K=l 



CONCLUSION 



1<K<2 



K=2 



We have studied in this paper a one-dimensional sys- 
tem of spinless fermions submitted to a quasiperiodic 
potential. Using a bosonization technique to treat the 
interactions exactly and the RG approach we had in- 
troduced in Ref. 12, we have investigated the effects of 



FIG. 9: A possible evolution of the density of states as a 
function of the Luttinger liquid parameter K. For K = 1 
the spectrum is a set of measure zero and gaps are present 
almost everywhere. For 1 < K <2 the small gaps close as K 
increases. For K = 2 even the largest gaps are closed since 
even a periodic potential is irrelevant above K = 2. 



Indeed the frequency dependence of the conductivity 
become for u < ujp — u/^, ^((jj) ^ oj^ (up to log correc- 
tions) as for a non interacting system. 

In the quasiperiodic case, the strong couphng fixed 
point is elusive so we can make only educated guesses. In 
the noninteracting case {K = 1) for a point of the spec- 
trum, the correlation functions decay as a power-law. It 
is thus unlikely that for the interacting case the strong 
coupling regime has a characteristic lengthscale in a sim- 
ilar way than the commensurate or disordered system. 
Thus the most likely possibility is that for the quasiperi- 
odic case the length ^ = ae^" separates, two power-law 
regimes with different exponents. Since for the nonin- 
teracting quasiperiodic case, contrarily to the disordered 
case, the correlation functions are still power-law, inter- 
actions are likely to still play a role even in the strong 
coupling regime. One can thus naively expect in that 
case that the exponent in transport and other correla- 
tion functions still depend on the interaction strength, 
albeit probably in a different way than in the weak cou- 
pling regime. This crossover is schematically shown on 
Fig.|. 

The spectrum can also be inferred. For the non inter- 
acting case K — 1 it consists in an infinite set of zero 
width bands. Clearly, the largest gaps must closed for 
iiT = 2, as it is the case for the pe riodic system. Since 
there is, as shown in Section III C 2| , a MIT at K = 1 for 
some filling fractions, the smallest gaps should close for 
K — 1 + e {e -^ 1). We thus naturally expect larger and 
larger gaps to gradually close as the interactions become 
more and more attractive and the Luttinger parameter 
moves from K = 1 to K ^ 2. Such an evolution of the 
spectrum as a function of K is depicted in Fig. pi It 
would be interesting to check this scenario by numerical 
investigations. 



various types of quasiperiodic potentials known as pre- 
cious potentials. We show that quasiperiodicity leads 
to a novel class of MIT as a function of the strength 
of the interactions, since for special filling factor the 
transition is pushed to the non interacting point (insu- 
lator for repulsive interactions and metallic for attrac- 
tive ones). We have determined the critical exponents 
and the associated lengthscales and showed the univer- 
sality of the results for all types of precious potentials. 
Our results a re . in n gQ od agreement with recent numerical 
investigationalilOtj. 

We have also analyzed the transport properties such 
as the conductivity and the conductance. These quanti- 
ties behave as power-laws with respect to the tempera- 
ture (or with respect to the size for the conductance at 
T — 0) with an interaction-dependent exponent. In the 
metallic regime, the RG flow converges towards a fixed 
point which allows to extract the full properties of the 
system. In the insulating regime, the system flows to 
a strong coupling regime at a lengthscale that we have 
determined. Above this lengthscale, it is necessary to 
analyze the physics with non pertubative method. For- 
tunately, one can still estimate qualitatively the behavior 
of some of the quantities. For instance, concerning the 
spectrum one expects, as the interaction becomes attrac- 
tive (or equivalently as the Luttinger parameter reached 
K > 1), that the smallest gaps start closing until K — 2 
where all gaps are closed. Similarly for the tempera- 
ture dependence of the conductivity one expects below 
the crossover scale Tc, a power-law dependence of the 
conductivity. It would be interesting to check the aboite 
proposals in-the numerical solutions, both for spinlessEj 
and spinfulEa systems. More generally, the extension of 
these investigations to other types of potentials such p^ 
the Prouet-Thue-Morse or the paper-folding potentialsc2l 
would be very useful since they also display a complex 
(dense) Fourier spectrum. Being able to tackle strong 
modulations would also allow for a cpaiparison with po- 
tentials such as the Harper potentialal2l which at the per- 
turbative level are similar to simple periodic ones. 

Several other questions are prompted by our study. 
The first one concerns the temperature dependence of 
the conductivity. The formula we derived assume that 
there are phase breaking processes, so that the temper- 
ature acts indeed as a cut-off in the RGE3. The validity 
of this assumption has bp2n recently explicitly proven for 
the Mott (periodic) caseEJ. In that case, the phase break- 
ing is provided by higher order periodicity (higher order 
umklapps). In the absence of such terms, the conductiv- 
ity would remain infinite. It would be interesting to carry 
on the same type of memory matrix approximation for 
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the quasiperiodic case. Second, it would be interesting to 
know if a correlated disorder could induce the same type 
of MIT as the one encountered here for the quasiperiodic 
system. Indeed, this type of disorder is susceptible to 
also produce a non trivial Fourier spectrum and thus to 
have a critical K that depends explicitely on kp- 

To conclude, we address the question of experimen- 
tal realizations to directly observe these effects. One 
could think about quasiperiodic chains in various de- 
vices. Since in one dimension there is a difpft equiva- 
lence between spinless fermions and bosonaljO it is pos- 
sible to investigate the physical properties of quasiperi- 
odic chain in Josephson junction arrayscS. The advan- 
tage of such systems, besides the excellent control that 
one can have on the potential, is that we can reach the 
attractive regime. One could also malize a quasiperi- 
odic chain using-gpantum dot arraysCj or patterning of 
a quantum wireCJ'EH^ . Finally, it would be interesting 
to check whether quasiperiodicity is relevant to describe 
systems such as DNA faL|j¥|hich there has been recent 
transport measurementsaoEj. 
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APPENDIX A: FOURIER TRANSFORM OF 
PRECIOUS MEAN POTENTIALS 

We consider a periodic one-dimensional chain deco- 
rated by a diagonal potential W whose amplitudes on 
each site can take two discrete values Wa = +A/2 or 
Wb = —A/2 according to the substitution rule : 



A -^ A^B, B ^ A. 



(Al) 



To compute the Fourier transform of M^, it is convenient 
to consider the Z*^ order approximant of the potential 
obtained by iterating (/— 1) times the rule ( Al). We have 
represented below the first approximants of the Fibonacci 
sequence (fc=l) : 



(A2) 



For the V"^ order approximant, W is thus a periodic po- 
tential with an elementary period ni = Fi containing 
si = Gi elements Wa and pi = G;_i elements Wb where 
(F„)„gN' and (G„)„gN are the precious mean sequences 
defined by : 



/ = 1 


B 


Z = 2 


A 


1 = 3 


AB 


I = 4 


ABA 


1 = 5 


ABAAB 


1 = 6 


ABAABABA 



70123456701 
30527416305 



BABABAABAAB 

70123456701 
30527416305 

FIG. 10: A possible correspondence between the potential 
and the length applied to the (f^ approximant of the Fi- 
bonacci sequence (ne = 8, se = 5, et n'^ = 5). The upper 
line gives the canonical numbering of the sites and the lower 
line, the conumbering. 



Fn+l ■ 
Gn+1 



kFr, 
-kOr, 



Fn- 
~ Gn 



In the quasiperiodic limit {I 



lonEa 



Vn > 1 (A4) 

Vn > 0. (A5) 

the ratio si/pi con- 
(fc 



verges toward the Pisot solutionE^ t = {k + Vfc^ + 4)/2 
of the equation x^ = kx + 1. For k = 1,2,3,4, r is 
known as the golden, silver, bronze and chocolate mean 
respectively. Since r is always irrationnal (Vfc G N*), the 
length of the period becomes infinite, and the sequence 
is quasiperiodic. 

The precious mean sequences can also be built by the 
Cut and Project algorithm starting from the usual Z-^ lat- 
tice and choosing, for the cut slope, a/ = si/pi- In this 
case, one obtains a periodic structure with two types of 
lengths La and Lb distributed according to (Al) but 
with a different origin than that given by (Al). To es- 
tablish a correspondence between the potential and the 
structure, one can, for example, affect Wa (resp. Wb) to 
a site if the adjacent left segment of this site has length 
La (resp. Lb) as displayed in Fig. 10. 

For this type of tiling (codimension 1), it is possi- 
ble to label each site according to their local environ- 
ment and thus-to, classify them with respect to their on- 
site potentiaE2l'L3^ . This procedure, known as the con- 
umbering scheme, allows to easily compute the Fourier 
transform of the potential, since it is simply expressed in 
terms of the so-called generating vector. 

To achieve this conumering in our case, we introduce, 
following Ref.E^, the parallel space vector in the Z^ lat- 
tice for the l*'^ order approximant A{' = (G'/,G'/_i), and 
we look for the generating vector g; determined by the 
condition det(A| , g/) = l.j^ince, the precious mean se- 
quences verify the relationEj : 



GiGi-2-Gf_, = {-iy, (A6) 

the generating vector is straightforwardly given by : 

gz = (-l)'(G,_i,Gz_2). (A7) 



Fi = F2 = 1 ; Go = 1, Gi = 



(A3) 



After the projection step, each site can then be indexed 
by its conumber j defined by rj = jn[ [ni] for j = [0, n— 1] 
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where rj denotes the canonical indexing (see Fig. 10), 
where nj — (— l)'(G;_i+G';_2) and where n^ = G1+G1-.1. 
These two numbers (n/ and nj) are in fact the lengths of 
the vectors A{' and g/ respectively measured in the Z^ 
lattice unit. Thus, all the sites whose conumber j e 
[0,si — 1] (resp. j e [si,ni — 1]) have a potential Wa 
(resp. Wb)- 

As a result, the Fourier transform of the potential is 
simply given by : 



We also introduce the right and left Fourier components 
of the fermions density operators : 



pAi) = J2''l,k+q^'',k, 



which satisfy bosonic commutation relations 



[pR{-q),PR{q')] = — V' 



27r' 



W q 



27rm\ 
n J 



^ 71 — 1 



,iqrj 



j=0 



[PLi-q),PL{q')] = -2^V' 
[PRiq),PL{q')] - 0. 



(B5) 



(B6) 

(B7) 
(B8) 



n-l 




^iqjn [n\ 



The commutation relations of these operators with the 
Hamiltonian Ht : 



nsm I 



(A8) 



[Ht,PRiq)] = vpqpBiq), (B9) 

[Ht,PLiq)] = -vpqpUq), (BIO) 



for integer values of tti g [1, n — 1] and : 

WiO) = X{s-p)/2n. (A9) 

Note that in (Rq) , we have omitted the index / for clarity. 



APPENDIX B: BOSONIZATION OF SPINLESS 
FERMIONS 

Let us first consider free fermions, i. e. with V = 
W = 0. The kinetic part of the Hamiltonian (|l|) is easily 
diagonalized via Fourier transform : 



Ht 



E^(^) 



4cfe, 



(Bl) 



where cj. = 1/VL ^ e*'''-'ri and e(fc) — ~2tcosk. 

If one is interested in the low-energy properties of 
the system, the only relevant states are those standing 
around the Fermi points ±kp. One can thus linearize 
the dispersion relation around these points and obtain 
an effective Hamiltonian : 

Ht = VF ^(fc - I^F)c'ji^k C-R,fc - {k + kF)c\j, CL,k, (B2) 



where vp — 2i sin kp is the Fermi velocity and where 
we have introduced the right R (respectively the left L) 
movers fermions with momentum close to +kF (respec- 
tively to —kp)- 

We now introduce the fermions fields : 



i'uix) 



VI? 






(B3) 



where v — R or L, so that the Hamiltonian (B2) writes 
in the continuum limit : 

Ht = ~ivp / dx ■4'Rix)dxi^Rix) - ■ipl{x)dooipL{x). (B4) 



explicitely show that piy{q) generate eigenstates of Ht 
with the energy vpq. This allows to write the kinetic 
energy as a bilinear operator in the bosonic fields : 



Ht 



TTVp 



Y^ [pR{q)pR{-q) + PL{q)pL{-q)] 



or in the real space : 

Ht = TTVp I dx [pr(x) 



PL{xf 



(Bll) 



(B12) 



where Piy{x) — ip}j{x)'ipi,{x) for v — R 01 L. 
We now introduce the fields 6 and 9 : 



</)(x) = 

e{x) = 



-j^lsr^ {p_R{q) + pdq)) 



-iqx 



q^O 



q 



^TT Y^ ipRJq) - PLJq)) -igx 

T. 2^ n ^ ' 



9^0 



together with their conjugate momenta : 



1 



n0(x) = -dJix) 



Y^^iPRiq)- Phiq))^ 

q^a 
Pr{x) - pl{x), 

1 



~'iqx 



Ilg{x) = -dx(t){x), 

n 
= ^Y.^PR{q)+PL{q))e-'^'^, 

= -Pr{x) ~ Pl{x), 
which obey the canonical commutation rules : 



[0(x),n4x)] 

[0{x),ne{x)] 



iS(x — x), 
iS{x — x'). 



(B13) 
(B14) 

(B15) 
(B16) 

(B17) 
(B18) 

(B19) 
(B20) 



(B21) 
(B22) 
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This finally allows to rewrite the kinetic energy : 



Ht 



VF 

27r 



dx [{-KlVf + {d^ 



(B23) 



with n = n^, which is the Haniiltonian of a one- 
dimensional elastic string. The fermions operators, and 
thus, all physical quantites (charge density wave. Cooper 
pairs density,...), can be easily expressed in terms of the 
(f) and n fields. The correspondence between the two sets 
of operators are given by the following relations : 






\r2nux 
1 

v27ra 



(B24) 
u\^-{'i>{x)+6{x))^ (B25) 



a being an ultraviolet cut-off and Q(x) being defined by : 
Q{x) — 7r/_^ n(a;')rfa;'. The operators Ur and Ul are 
anticommuting operators which ensure the correct com- 
mutation relations of the fermions operators. Note that 



these operators give essentially vanishing contributions 
in the thermodynamical limit and can be safely ignored, 
at least in the case of single chain systems. 

The great advantage of this bosonic formulation is that 
the interaction term Hy can be almost trivially taken 
into account. Let us consider, for instance, a two-body 
interaction : 



^v = 1^ I dx 



dx'p{x) V{x ~ x') p{x'), (B26) 



where p{x) is the fermion density operator on site x. 

If one is only interested in the low-energy physics, the 
fermion field operator il}{x) writes : 



^{x) 



ikpx 



'>Pr{x) 



— ikjpx 



i^dx), 



(B27) 



which correctly describes the low-energy behavior of the 
particle-hole and particle-particle density state. It fol- 
lows that the density operator expresses as : 



p[x) = i)\x)il}{x), 

= i^R{x)ipR{x) + ■iIjI{x)iPl{x) + e 

= PRix) + pl{x) + e-2^'=-^^l,(x)V'L(x) + e2^^-^^i(x)Vi?(x) 



-2'"'--V1j(^)^l(^) + e2^'=--Vl(2:)V'fl(x), 



-—dx<p{x) -\ cos(20(x) — 2kpx), 

n 7ra 



(B28) 
(B29) 
(B30) 

(B31) 



r 



and : 



p{x)p{x') = Pr{x)pr{x') + pl{x)pl{x') + pr{x)pl{x') + pl{x)pr{x') 



-i2kfix-x') /t 



i^],{x)MxWL{x')^R{x') + 



J2kF(x-x') /t I 



e— - - ' i!l{x)il}R{x)il}\{x')ijL{x') + 
oscillating terms. 



(B32) 



r 



In the case of short-range interactions the remaining 
oscillating terms are supposed to give negligible contribu- 
tions and we keep only terms behaving like e"'^*^^*^^"^ ). 

If one consider local interactions, one can write the 
interaction term in the general form : 



i/v 



jdx 



Pr{x)pr{x) + pl{x)pl{x) 



i± 
2 

92Pr{x)pl{x) 



(B33) 



with the standard notations of g—ology. Gathering (B23|) 



and (p333| ) and the Hamiltonian HQ{t, V) then becomes : 

Ho{t,V) = ^Jdx {uK){7TUr + (-|) (a,0)2, (B34) 

where we have introduced the parameters u and K (Lut- 
tinger parameter) given by : 



K 



'(•-f)' -(!!)• 


(B35) 


/ 2'kvf + .g4 - .92 
2ttvf +94 + 92' 


(B36) 
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The Hamiltonian (B34) can still be interpreted as the 
Hamiltonian of an elastic string with effective parameters 
u and K. 

For the model considered here, the interaction term 
writes : 



Vy^niUi+i -^ limVa j dx p{x) p{x + a) . (B37) 



the interactions are strong provided the correct u and K 
parameters are used. 



For example for kp — 7r/2 one has for the interaction 
( p37| ) : 



This expression can also be written as (B33) for 



F < i with : 54 = 32 = 2V^(1 - cos(2fcFa)). Note 
that, in the non interacting case one has K = \ and 
u = vp — 2isin(fci?a). For repulsive interactions 
{V > Q), K < 1 while the attractive case (T^ < 0) leads 
t o K > 1. A remarkable fact is that. t|hp. representation 
( |B34|) is, in fact, completely generaOCJ and gives the 
correct low-energy description of the system, even when 



V_ 
2i 



2Kt . f f 1 

sm TT 1 

2K -I \ \ 2K 



(B38) 
(B39) 



Concerning the bosonization of Hw, one straightfor- 
wardly obtains, in the continuum limit : 



H 



w 



dxW{x)p{x) = dxW{x) ( dx(f>{x)~\ cos(20(a;) 



2kpx' 



(B40) 



r 



APPENDIX C: RENORMALIZATION GROUP We evaluate the correlation fimction R{x,t,x' ,t') per- 

EQUATIONS turbatively at second order in power of the coupling con- 

stant g : 
The aim of this appendix is to derive the RG equations 
used in Section II, following the method of Ref. 54 



J 



TrC 



i^/2<p{x' .t') —i^/2(f)i^x^T) 



= h- 



9 



u(27ra)2 
I 



/■ + i 



2 u2(27ra)4 



0(5^) 



(CI) 



with 



Tt-C 



i\/2(j){x' ,r') —i\/2(p{x.r) 



h ^ J2 ^(Q) / '^'^i (^- 



^i\/2<p(x' ,t' ) -i\/2<p{x,T) 



g4[20(xi,ri)+Q xi\ _|_gi[-20(xi,ri)+Q + xiJ 



with Q* ~ Q ± 2kF and 



(C2) 
(C3) 



72= ^ t^(Qi)Ty(Q2) /"d^n /"d2r2/r^e*^'^(^''^')e-'^^(^'^) X (C4) 

l^i[24,(xi,Ti)+Q-xi] _|_ gj[-20(i:i,Ti)+Q+a:i]l I ^i[2^{x2,T2)+Q- X2] _|_ ^i[-24,{x2,T2)+Q+ X2]'] ^ 
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where ()o denotes the average performed with respect to 
the free action, i. e. with g = 0. The index c in Eq.(C4) 
refers to the connected correlation function with respect 



J 



to the quantities indexed by 1 and 2 



The average in Eqs. ( |C2[ ) ( |C3| ) and ( |C4[ ) are easily 
performed using the relation : 







I In ■ 



(C5) 



r 



when \ri — rj\ ^ a. Note also that the expression ( |C5| ) with F{r — r') = Kin |r — r'|/a. 
vanishes if the "neutrality" condition, ^^ c,; = 0, is not 
satisfied. One thus has : 



/o = (Tre 



^iV2(j>{x' .t') —iV2il>(x,T) 










(C6) 
(C8) 



Note that Ji vanishes since it does not respect the 
neutrality condition. It follows that the first non triv- 
ial contribution of t he potential is provided by the term 
/2. Keeping in (C4) the terms satisfying the neutrality 
condition and performing the average over the fields one 
obtains : 



E W{QimQ2) I d?v, I d'r, e-2^(---)g^^ /:; Wiy).y 



(C9) 



-V2F(r-ri) + V2_F(r-r2)+V2-F(r'-ri)-V^F(r'-r2) _ ^ 



pjiQl Xi+Q^X2) j^ pi{Qtx2+Q2 xi) 



At this order the potential enters only through a pure 
phase and does not drastically affects the long distance 
behavior. Anyway, a choice of a vanishing average value 
for the potential renders this effect negligible. 

The integration over ri et r2 is performed with the 
following variable changes : 

St St 

ri=R+y, r2 = R-y, (CIO) 



since, due to the factor e ^^(''i ''2) in Eq.(|C9|) the non 
vanishing contributions are given by "points" separated 



by small values of Sr = Vi — r2. If the potential is trans- 

ilK. pi W(y)dy 

lation invariant, the term e " ^2 does only de- 

pends on xi — X2 and is thus of order one for small 5v. Of 
course, in the quasiperiodic case, this is no longer true. 
Nevertheless, if one is only interested by the intrinsic 
properties of the potential, it is reasonable to perform 
an average over all the possible choice of the sequence 

i2ii Pi W{y)dy 

origin. Then, we can also consider that e " -'^2 
only depends on xi — X2 and replace it by 1. 

At leading order in 5r, one has : 



F{v - ri) - F(r - V2) ~ F{v' - ri) + F(r' - rs) = <5r.VR [F(r' - R) - F(r - R)] + o{5rf . (Cll) 

I 



On the other hand, with R = {X,Y) et 5t — {6x,Sy) and 
one has : 



Sx 



ox 
Q^xi + Qtx2 = {Qi +Q2)X+ {Ql - gj) — , (C12) 



}tx2 + Q2X1 = {Qi + Q2)X + {Q2 - Qt) ^, (C13) 
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which imphes that the integration over the variable X Qi and Q2- Also, taking the fact that the potential W 
selects the contributions Qi = —Q2 = Q in the sum over can be taken as real one has W(Q) — W*{—Q) and : 



I, = e'^/.X^)^^e-^(-^-'-')5]|w^(Q)|Yrf'R/'^''5r e-^^^^r) 



j{Q-Sx) _^^~~i{Q + 5x) 



X (C14) 



<5r.VR [F{r' - R) - F(r - R)] 



r 



Let us now define : C — F{^' ^ f^) ^ ^i^ ~ R-)- 0^.6 has : 

(Jr.VR 0' = {Sx dx 0' + {Sy dy C? + ^Sx 6y dx C dy C- 

(C15) 
By parity only the two first terms would survive to the 
integral over 5y (or Sx). On the other hand, the integral 
over R is easily performed and leads to two kind of terms : 



fact that the free correlation function F writes, in fact, 
as : 



F(r- 



Kln^- — ^+dcos(26lr-r'), (C18) 



d^R C {djc + d^)C^ -AnK^ In 



r — r 



(C16) 



and 



f dRC{d'^x-dl)C = -27rcos29r-r', (C17) 

where ^r-r' is the angle between the vector (x, r) and the 
X axis. The occurence of the second term comes from the 



where d parametrizes the anisotropy between the space 
and time direction. One ha s d — in the original Hamil- 
tonian but, as seen in Eq.( C17| ), this anisotropy is gen- 
erated by renormalization. Ultimately, this is equivalent 
to a renormalizaton of the Fermi velocity u. However, 
since the anisotropy itself is of order g^ , the correction to 
u can be neglected in this calculation. 

Gathering the preceeding results one thus finds : 



J 



I, ^ e'"^/>(^)'^^e-^("-'-')2^ir2^|w^(Q) 



2 \v-v'\ 
In-i X 



(C19) 



d5r<52re-2^(*'-) 



^i{Q-&x) , ~i{Q-^&x) 



which also reads 



I2 ^ e'^/.>(^)^^e'^("^-'-')2^ir2^|w^(Q) 



Q 



In 



r — r 



+00 
lixoi^ I ddr 



3-2Kr 



MQ^Sr)+JoiQ-Sr) 
I 



(C20) 



where Jq is a Bessel function. 



Finally the correlation function R{x, r, x' , r') writes : 
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/j,^^iV2ci>(x'y)^~i^,j,(x,T)\ ^ g»^ /J; W{y)dy g-F(r-r') K ^ 



^2^2 |j. _ J., 



In^^^kl^(Q) 



STT^Tt^ a 



(C21) 






3-2K 



Jo(Q"fr)+Jo(g+<5r) 



By re-exponentiation one has 



R{x, r, a;', r') = 



ATe 



(C22) 



with 



^e//=if-^5]jy| 



r+°° d5r 



3-2K 



3Q{Q+6r) + 3o{Q-Sr) 



(C23) 



with fjQ — \aW{Q)/u. 

To derive the RG equation, one has to consider an in- 
finitesimal variation of the running cut-off a{l) to a{l) 



,dl 



in Eq.(C23). This leads to 






with 



s(0 = lT.yl [Jo(Q+«(0) + Jo(g-«(0)] 



(C24) 
(C25) 

(C26) 



Here one has to note that the use of a sharp cut-off in real 
space leads to the occurrence of Bessel function Jq . This 
choice is in fact not satisf ying s ince it does not ensure the 
convergence of the sum ( |C26| ) . This is the reason why 
we have considered more general cut-off procedures for 
which Jo is replaced by faster decreasing functions which 
typically satisfy (\ 



;ed by 

HO) 



APPENDIX D: MEMORY FUNCTION 

In a normal metal (finite conductivity at w = 0) the 
Kubo formula : 



a{uj) = - 



2uK 



+ x(^) 



(Dl) 



implies that x(0) = —2uK/Tr. Then (Dl) can be reex- 
pressed in terms of a meromorphic function M through : 



a{uj) 



i2uK 



1 



IT CJ + M{ljj) ' 



(D2) 



where M is given by 



M{uj) 



ujx{i^) 



X(0) - x{^) 



(D3) 



The interest of the function M lies in the fact that, con- 
trarily to the conductivity itself, one can expect M to 
have a well behaved expansion in the scattering poten- 
tial Hw- Indeed in a simple hydrodynamic approxima- 
tion M{ijj — > 0) would simply be the inverse relaxation 
time M ^ i/r, leading to the standard Lorentzian broad- 
ening of the Drude peak. Another way to formulate it is 
that a perturbative calculation of the memory function 
is close to a perturbative calculation of the resistivity. In 
the lowest order in the scattering potential Hw one gets : 



and : 



X(0) - xi^) ^ X(0), (D4) 



u;x{u;) = [{F;F)'^-{F;F)l^,]/u;, (D5) 



where F operator takes into account the fact that the cur- 
rent is not a conserved quantity F — [j, H] and {F; F)^ 
stands for the retarded correlation function of the op- 
erator F at frequency lo. Since F is itself proportional 
to the scattering potential, at lowest order the average 
can be computed with the Hamiltonian in the absence of 
scattering potential. This leads to : 



M{uj) = 



[{F;FZ-{F;FZ^,]/u; 

-x(o) 



(D6) 



where O*^ stands for an average with (pi) only. Since all 
averages are to be computed with the quadratic Hamil- 
tonian (H) only, the computation of M is now feasible. 
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